Dispersion interaction between crossed conducting wires 
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We compute the T = OK Van der Waals (nonretarded Casimir) interaction energy E between 
two infinitely long, crossed conducting wires separated by a minimum distance D much greater than 
their radius. We find that, up to a logarithmic correction factor, E oc —D~^ Isin^p'^ f{9) where f{9) 
is a smooth bounded function of the angle 6 between the wires. We recover a conventional result 
of the form E oc — Z)"'' \sm9\~^ g{9) when we include an electronic energy gap in our calculation. 
Our prediction of gap-dependent energetics may be observable experimentally for carbon nanotubes, 
either via AFM detection of the vdW force or torque, or indirectly via observation of mechanical 
oscillations. This shows that strictly parallel wires, as assumed in previous predictions, are not 
needed to see a novel effect of this type. 
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At the micro- and nano-scale, dispersion (van der 
Waals, vdW) forces are ubiquitous [1], and recent ad- 
vances in manufacturing and measurement techniques 
have prompted much interest in their precise form. 

The simplest theories sum vdW interaction energies 
between pairs of molecules, which is a good approxima- 
tion for dilute insulating objects, where the dipole fluc- 
tuations at different points of one body are almost inde- 
pendent. For non-dilute dielectric and magnetic mate- 
rials this summation approximation can be misleading, 
and may even give the wrong sign of the interaction 2]. 
Moreover, for anisotropic conducting nanostructures, the 
non-locality of Coulomb screening and associated density 
correlations within each object may change the form of 
the dispersive forces altogether 0, |3| . 

This physics is exemplified by the class of quasi- ID ob- 
jects, which exhibit correlation phenomena of both theo- 
retical and experimental interest. Indeed, in the extreme 
limit truly confined ID electrons experience a Luttinger 
liquid instability, as (e.g.) in single walled armchair nan- 
otubes 0. 

When two such "wires" are placed parallel and close to 
each other, the coulomb interaction between their density 
fluctuations may become a relevant perturbation, result- 
ing in rich behavior at low temperatures and densities, 
such as locked charge density waves and Wigner cristal- 
lization. Q. The density density interaction is also re- 
sponsible for coulomb drag phenomena whereby a cur- 
rent applied to one wire induces voltage on the other 
wire I 

Dispersion forces between ID systens are interesting 
even at larger separations. For the case of two in- 
flnitely long strictly parallel wires separated by distance 
D greatly exceedingthcir radius h, the interaction energy 
is known [3, lloL llll. Q| to be strongly dependent on 
the presence of an electronic energy gap: 



£""^^11 cx -D-'^{\n{D/h))-^''^ {metallic). (1) 
E'"'^^/L oc -D-'^ (semiconducting) (2) 

The result (P) was obtained via zero-point energies of 
coupled plasmons in the random phase approximation 
(RPA), followed by a perturbative evaluation of the re- 
sulting integral Li. The form d) was also supported 
by diffusion Monte Carlo calculations • 

Such forces may be important when one considers so- 
lutions of nanotubes or long molecules. Indeed, while 
"wires" minimize their energy by aligning, in a solution 
of nanotubes a parallel configuration might not be formed 
because of entropic reasons. Therefore, it is important to 
understand the wire-wire interaction for a general orien- 
tation. Here we consider the angle and distance depen- 
dence for wires that are well separated. We will show 
that the (absolute) vdW energy of a pair of non-parallel 
wires inclined at angle 9 is, up to a logarithmic correction 
factor specified later, 
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f{e) (metallic) (3) 
g{9) (semiconducting) (4) 



Here D is the least distance between points on the two 
wires, and 9 is the angle between the wires. / and g are 
smooth bounded functions. In this limit there is also a 
prospect of measuring the dispersion force between two 
nanotubes directly, via Atomic Force Microscopy (AFM) 
or spectroscopy of mechanical vibrations. Eqs p|4p show 
that unusual gap-dependent results are expected from 
such experiments not only when the tubes are parallel 
but also when they are non-parallel. 
The crossed-wire interaction has previously been re- 
lated to the interaction between anisotropic media in 
a similar manner to the way Casimir forces between 
molecules can be obtained from the Lifshitz formula by 
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taking the dilute limit. The interaction between non- 
isotropic materials related to our problem was studied in 
[ist : there, the energy of a pair of media conducting only 
in p rescribed (but different) directions is computed. In 
|14| the interaction between non-isotropic dielectric me- 
dia was considered and the limit of a dilute medium was 
associated with the interaction between a pairs of wires. 
The asymptotic results of IJ] were similar to those of a 
pair summation approach (i.e. like Eq. (jl])) even for 
metallic cases), and qualitatively different from Eq. ^) 
derived below for metallic wires. We suspect that this 
difference may be due to the different Coulomb screen- 
ing physics of a single pair of metallic wires compared to 
an infinite array of such wires. (See [l6| for similar con- 
siderations relating to layered systems). This question 
is related to the very non-additive dispersion physics of 
low-dimensional, zero-gap systems in general (Sjj . 

In 17[, the orientational interaction between 
spheroidal dielectrics has been computed, but the 
metallic limit treated here was not yet considered. Some 
related work has also been done on ID conductors 
in a coUinear "pointing" configuration [l^, but the 
calculation appropriate for metallic ID conductors does 
not appear to have been done for cases where they are 
aligned in a non-parallel geometry. 

Here we will compute the dispersion energy directly 
for crossed wires within second order perturbative non- 
retarded theory, as justified below. The interaction en- 
ergy in this approach is (see e.g. [l^) 
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~p',iu)x^^\p,p',tu) (6) 



Here x(r, r') is the density-density response of each sub- 
system in the absence of the other sub-system, and 
x(PiP') is its 3D Fourier transform. V{r) is the bare 
Coulomb interaction between electrons in the two subsys- 
tems, and V{p) = 47re^/|p|^ is its 3D Fourier transform. 
V can be considered small in the present case because of 
the large inter-wire separation assumed here. Indeed it 
is easily verified that (O yields H]) for the case of parallel 
well-separated ID conductors. 

We consider a pair of wires A, B described in FiglT] 
One wire is assumed to lie along the z axis. The xz 
plane and the origin are determined by demanding that 
the point on wire B, lying closest to wire A, is the point 
(i?, 0, 0). The orientation of wire B is then determined by 
a single angle 9. 

For later convenience we introduce unit vectors in the 
yz plane, parallel and perpendicular to wire B: u = 
(O,sin0,cos0), V — (O,cos0, — sin0). Locations along 
wire A are denoted z, and locations along wire B are 
determined by a signed ID position variable s so that 




FIG. 1: Van der Waals interaction of two crossed wires A and 
B at distance D and angle 9 



general points on wires A and B are TIa{z) — (0,0, z), 
R-s('S) = Dx + su 

The density-density response functions in 3D space of 
the two wires are written in terms of an assumed strictly 
ID density-density response function x(2;,a;) for electron 
motion along a wire: 

x'-^\v,v',oj) =x{z-z',u)5{x)5{x')5{y)5{y') (7) 
X*^ni",i-',w) = x((i" - r') • u,w) 

x5{x - D)5{x' - D)5{y ■ v)5{v' ■ v) (8) 

We express x(-Z)'^) in terms of its ID Fourier transform. 
xiz^oj) = ^ J x(9j '^) exp(igz)(iz. Then, doubly Fourier- 
transforming ([7]) and ([5]) we obtain 



x(^)(q,q') =27r/dfcix(fci,c.)%, + fci)<5(<zi-fci) , (9) 
X*""^ (q, q') = 2^ / dk2x{k2,u;) exp(i(fe + q'jD) 

xS{cos9qz — siiiOqy + k2)S{cosOq'z — smOq'y — 

(10) 

Putting ([grrUI into dS]), and defining 

F{ki,k2) = / x{ki,iu)x{k2,iu)du (11) 



(12) 



we obtain 
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X exp(i(pa; +p'^)D)5{ki -pz)5{p'z + ki) 
xS{cos6pz — sinOpy + k2)S{cos9p'z — sinOpy - 
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(13) 



We calculate the single-wire response x for eqs (|11I13[) 
within the Random Phase Approximation. From the 
work of Li and Das Sarma [20|, this RPA approach 
for the long-wavelength collective motions can even be 
applied to a strongly interacting model of a ID con- 
ducting system such as a Luttinger liquid. This would 
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not apply for broken-symmetry cases such as a Wigner- 
crystalline wire, but this case is likely to occur only for 
small inter- wire separations which are not the focus 
in the present work. 

The starting point of the RPA for a single wire is 
the bare density-density response xoik^uj). The simplest 
model for xo j at small k and small lu = iu, is 



Xo(fc,«u) 



noiD 



(14) 



where noio is the number of electrons per unit length 
in the ground state, and mis an appropriate mass, loq 
(with hwo w Bloch energy gap Eg) is a harmonic pinning 
term allowing the modelling of both semiconducting and 
metallic wires. We take ujq ~ for the metallic case, and 
ujQ for semiconductors. (|14p follows from classical 
motion of individual independent electrons (F = ma) 
but is also valid for the quantal motion of independent 
Fermions at small wavenumber fc, corresponding, in the 
present case, to large interwire separation D ~ k~^. 

The Coulomb interaction e^/|zi— ^2! between elec- 
trons on a strictly one-dimensional wire has a divergent 
Fourier transform corresponding to the region zi —>■ Z2. 
In real quasi-lD structures such as carbon nanotubes 
there is a finite radius b corresponding to the lateral 
spatial extent of the one-electron orbitals. This serves 
to smear the Coulomb divergence at zi Z2. The re- 
sult depends on the detailed cross section assumed. For 
definiteness we choose a smeared Coulomb potential for 
electrons free to move on a cylindrical shell of zero thick- 
ness and radius b, representing a nanotube. The smeared 
intra-tube potential for this case is analytic: 
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2e^\ln{\kb\)\, for\kb\ « 1. 
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where Iq, Kq are modified Bessels. The long- 
wavelength form (|16p is universal to all wire profiles but 
causes spurious collective modes and a divergent energy 
at larger k, requiring a cutoff to be imposed. Although 
the final result is highly insensitive to this cutoff, we pre- 
ferred to use the specific form (fTSl) . thereby avoiding a 
cutoff altogether. 

We now use the RPA to write the long-wavelength 
response of mutually interacting quasi-lD electrons to 
an external potential as 
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Here cji£i(fc) is a quasi-acoustic ID plasmon frequency, 
which can be written in terms of a ID velocity cuj as 

i^ioik) = ciD \k\ ^/lo{\bk\)Ko{\bk\), cid = \/2e^nmDlm 

(19) 



Now the ID plasmon group velocity dum / dq from (fT9|) 
is much less than the speed of light c except for q < 
0{b^^exp{—c^/cl^)). Thus electromagnetic retardation 
can be ignored except at extremely large distances D » 
6exp(cVc2^). 

Puting dni) into find 
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Case 1, conducting wires, — 0. 
To analyze ()13|) we transform to dimensionless plane 
polar coordinates {K, (f) : 

f) 

ki = D-^Ksin{(t)+ -), k2 = D-'^Ksin{(j)- -) (21) 

with Jacobian J = D^'^K/\sin{0)\. This yields 
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h{e,r), r ^ b/D » 1 (22) 
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with s± = |sin((/) ± |)| and i± = \J Iq {Krs±) Kq {Krs± ) . 
Note that h{6^ r) is well behaved as ^ 0, with a limiting 
value /i(0,r) = {ln{D /b))'^''^ for r = b/D «< 1. 
For r <<< 1 wc thus obtain the analytic result 

^'''"-l^^^^^M^( (ln(i^)5))3/0 ""'^^^ 

(24) 

For arbitrary 6 and r = b/D << 1 the function h{9,r) 
defined by (j23[) requires numerical investigation. We de- 
fine a smooth function A with correct angular period: 

N 

A{e,r) = h{B,r)(ln{\/r)fl'^ « ^ c„(r)cos(2n61) 

r!=0 

Cn{r)^ ^^f'f . , x = l/lnr (25) 

The choices iV = 5,/ = 3,J==2 give a good fit to ^ 
for X < 0.4, i.e. for D/b = r^^ > 16.5, and the corre- 
sponding coefficients a, b are given in Table 1. 

Case 2, semiconducting wires, luq ^ 0. Here there are 
two analytic cases according to the separation D: either 
of the terms on the right of could dominate for the k 
values k « D~^ that dominate the energy integral 

Case 2a, smaller separations When 



D «Do = ciij^o ^y/ln(ciD/ii^ob)) 



(26) 



we can again ignore ujo in (I18|) . thus recovering the 
"metalhc" vdW energy (f22|) . Because of our large- 
separation approximations, this conclusion only holds 
provided that Dq from satisfies Dq >> b, which 
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din ri=0 n=l n=2 n=3 n=4 n=5 

i=0 0.85708 0.11770 0.01534 0.00453 0.00189 0.00097 

i=l 0.18049 0.19040 0.04743 0.01538 0.04505 0.00228 

i=2 2.85682 0.51752 0.12230 0.03615 0.01378 0.00783 

i=3 2.78324 0.61964 0.13410 0.03497 0.03237 0.01795 

bjn n=0 n=l n=2 n=3 n=4 n=5 
j=l -2.29966 -0.24097 -0.06158 -0.00820 -0.00657 -0.00357 
j=2 2.47243 -0.19977 0.70113 0.91201 0.95313 0.97461 

TABLE I: Coefficients for log-cosine expansion (|25|) of A(^, r) 
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FIG. 2: vdW Force contributions in picoNewtons between 
(5,5) CNTs (1st 4 columns) and Coulomb force between two 
electrons. 



can occur for very small gaps loq such as that noted for 
(9,3) nanotubes in 

Case 2b, larger separations D » Dq Here we can 
ignore the ojpj^jj term in (fTS]) for the relevant values k w 
1/D. Then HID becomes F « 
(PTjl . we evaluate (fTS)) as 



^1^2?^ and using 



j^{2):seraicond 



Stt 1 



hc\jj 2 cos^ ( 
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(27) 



This is consistent with Eq (57) of [14 1. 

Conclusions: Eq (|22|) is the principal result of the 
present work, along with ([23]) and (|25p. It shows that 
nonparallel conducting wires experience a vdW attrac- 
tive energy that decays much more slowly with distance 
D than the standard D~'^ dependence predicted by sum- 
ming contributions over all elements i,j of the wires. 
(|22p also shows a strong angular dependence, giving rise 
to significant slowly-decaying vdW torques. 

The interaction predicted here might be measured di- 
rectly in Atomic Force Microscopy experiments on metal- 
lic carbon nanotubes, or indirectly via their mechanical 
vibrations. In Figure 2 we estimate the force F"^'^^ be- 
tween the conduction-band electrons of two freestanding 
metallic (5,5) carbon nanotubes in vacuo, as a function 
of separation Z), both for strictly parallel tubes of length 
1 micron using Eq. (1) of and for infinitely long 
tubes at angle 6* = 1° from Eq. F"'" oc D'^ is 

the dispersion force from the remaining insulating elec- 
tron bands via the usual pairwise summation approach 
using data from 



2l| and [22] ■ For comparison, F'^ 



gives the force between two localized (impurity) single- 
electron charges separated by D. The dispersion force 
pdzsp ^ pmet _^ ptns ^^^j^ distinguished from that 

due to any localised extra electrons because F'^'^^p, un- 
like F'^~'^, will be invariant when one tube is made to 
slide along its own length. 

We note finally that the present theoretical results can 
be understood in terms of long-wavelength collective ex- 
citations, and are not limited to T = QK. However in 
practice T needs to be low enough for a long electronic 
mean free path > Z? to be maintained. 
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